A = [3,-1,-1;-12,0,5;4,-2,-1];
x = [1,1,1]';

[eigx, eigA] = eig(A)

for i = 2:10
    x(1:3, i) = A * x(1:3, i - 1);
end

x